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Abstract 

The distribution of the initial short-time displacements of particles is 
considered for a class of classical systems under rather general condi- 
tions on the dynamics and with Gaussian initial velocity distributions, 
while the positions could have an arbitrary distribution. This class 
of systems contains canonical equilibrium of a Hamiltonian system as 
a special case. We prove that for this class of systems the nth order 
cumulants of the initial short-time displacements behave as the 2n- 
th power of time for all n > 2, rather than exhibiting an nth power 
scaling. This has direct applications to the initial short-time behav- 
ior of the Van Hove self-correlation function, to its non-equilibrium 
generalizations the Green's functions for mass transport, and to the 
non-Gaussian parameters used in supercooled liquids and glasses. 
KEY WORDS: Particle diffusion, non-Gaussian effects, cumulants, 
time expansion, Van Hove self-correlation function, Green's functions, 
supercooled liquids 
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On the structure of this paper 

This paper concerns a universal property of correlations of the initial short- 
time behavior of the displacements of particles for a class of ensembles 
of classical systems, both in and out of equilibrium (the latter one be- 
ing somewhat restricted). Among other things, these correlations (ex- 
pressed in terms of the so-called cumulants) have applications to neu- 
tron scattering^, 0, 0, 0, 0, 0, 0, 0, E3, E3, Q, to the description 
of a restricted class of non-equilibrium systems on small time and length 



scales 0, Q, 0, [Til. 17 , and to heterogeneous dynamics in supercooled 
liquids and glasses [[1, El, E3, EE . 

The property can be formulated as a mathematical theorem, which may 
have more applications than are considered in this paper. Our interest is 
however mainly in its physical applications. Since the general theorem and 
its proof are best formulated in a rather abstract way which is not necessary 
for the currently known physical applications, the paper is split up into a 
physical part (Part ^) containing the physical formulation of the Theorem 
and its applications, and a mathematical part (Part|2J) containing the general 
mathematical formulation of the Theorem and its proof. 

These two parts are further subdivided into sections. In Sec. 11.11 the 
physical motivation for studying particle displacements is given. In Sec. 11.21 
we introduce the physical systems that will be considered. Section H~3l gives 
the necessary definitions to be able to treat the cumulants of the particle 
displacements. We state the Theorem in physical terms in Sec. 11.41 and 
discuss its physical applications in Sec. 11.51 

As for the mathematical Part [21 of the paper, in Sec. 12.11 we give the 
general definition of cumulants of any number of random variables and dis- 
cuss some of their properties, while in Sec. 12.21 we present the Theorem in 
its full mathematical generality. In Sec. 12. HI we prove the Theorem. In the 
proof we need an auxiliary theorem concerning Gaussian distributed vari- 
ables whose proof is postponed to the Appendix. Some coefficients occurring 
in the Theorem are worked out in Sec. 12.41 

We end with conclusions, which are followed by the Appendix. 



1 Physical Part 
1.1 Introduction 

Our motivation to consider the displacements of particles originates from 
trying to describe the behavior of non-equilibrium systems on all time scales 
using the so-called Green's function theory, introduced by J. M. Kincaid[lj]. 
This theory aims to describe, among other things, the time evolution of 
the number densities, momentum density and energy density, by expressing 
them in terms of Green's functions. It has so far been successfully applied 
to self-diffusion [lH] and to heat transp ort 113 . while a study of mass 

transport in binary (isotopic) mixtures [17| is in progress. The main advan- 
tage of the Green's functions over hydrodynamics is that one expects that 
they can describe the system on more than just long time and length scales, 
in particular also on time scales of the order of picoseconds and on length 
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scales of the order of nanometers. The connection between the picosecond 
and nanometer scale can be understood by realizing that with typical veloc- 
ities of 500 m/s, a particle in a fluid at room temperature moves about 0.5 
nm in 1 ps. Hence studying the picosecond and sub-picosecond time scales 
could also be relevant for nanotechnology. 

For mass transport in multi-component fluids 3, 3h the- Green's func- 



tions G\(r,r',t) have the physical interpretation of being the probability 
that a single particle of component A is at a position r at time t given that 
it was at a position r' at time zero. At time zero the system is not in equi- 
librium, in fact, a class of far-from-equilibrium situations has been studied 
in this context [bll. Hfl HfiL Ht1| . Given this interpretation, it is clear that the 
displacements of single particles are the central quantities in the Green's 
function theory. 

Apart from the non-equilibrium aspect, the above interpretation of the 
Green's functions is the same as that of the classical equilibrium Van Hove 
self-correlation function G s (r — r',t)Q, U In the literature on the clas- 
sical Van Hove self-correlation function in the context of neutron scattering 
on an equilibrium fluid 0,0,0], the so-called cumulants of the particle dis- 
placements (defined in Sec. EH below) have been studied by Schofield[]| and 
Sears 0], among others. The relevance of the cumulants for the Van Hove 
self-correlation function can be seen from its Fourier transform, the incoher- 
ent intermediate scattering function [H. Ifj. IfiL Hnl| . which can be measured and 
is defined as 

F s (k,t) = (gifek-In^-nCOlX /ifcAx l( M ^ (L1) 



/ \ / 
/ eq \ / eq 

where k = /ck is a wave vector with length k and direction k (a unit vec- 
tor), ri(t) the position of a (single) particle at time t and () eq an equilib- 
rium average. In the last equality we have chosen k = x (in an isotropic 
equilibrium fluid the result is independent of the direction k) and de- 
fined Axi(t) as the displacement of the single particle in the x direction: 
Axi(t) = x • [n(t) - ri(0)]. In probability theory [^J Ei] log(exp[ifcA]) is 
the cumulant generating function for the random variable A, so we see from 
Eq. (jl.ll) that log F s (k,t) is the cumulant generating function of Axi(t). 
Note that Axi(t) is a random variable here as it depends on an initial 
phase point drawn from a probability distribution (here the equilibrium dis- 
tribution). The cumulant generating function is, b y d efinition, equal to 
YlnLi K n{ik) n /n\ where K n is the nth cumulant |23|, |2J]. The relation of 
the incoherent intermediate scattering function and the cumulants of the 
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displacement is thus expressed by0, 0, 0, E 



oo 

F s (k,t)=ex V ^2^(ik) n . (1.2) 

n=l 

This connection with the incoherent scattering function (and thus with neu- 
tron scattering 0, 0, 0, 0, 0, 0, 0, 0, 0, E3, E3, Q ) explains the early interest 
in the cumulants of displacements from a physical perspective. We note that 
Eq. (jl.2|) also shows that the cumulants could be obtained experimentally by 
measuring F s (k,t) via small k ("small angle") incoherent neutron scattering 
and fitting the logarithm of the result to a power series in k. In Sec. 11.51 we 
will discuss also more recent physical applications of the cumulants such as 
the non-equilibrium Green's functions and non-Gaussian parameters used 
in the theory of supercooled liquids and glasses. 

In the context of incoherent neutron scattering, it has been noted by 
Schofieldl and SearsQ that the cumulants n n of the displacement of a 
particle in a time i in a fluid in equilibrium with a smooth interparticle 
potential behave for small t as K2 = 0(t 2 ), K4 = 0(t 8 ), kq = 0(t 12 ), while 
all odd cumulants vanish. These results suggested for equilibrium systems 
with smooth potentials a behavior as the 2ra-th power in t for K n when 
n > 2 in general, but to the best of our knowledge no proof of this property 
is available at present. 

We note that an 0(t 2n ) behavior would be in stark contrast to results 
obtained for hard disk and hard sphere fluids in equilibrium [?1 E^l] . Sears 
considered results up to the 8-th order in t and the 12-th order in t for K4 and 
kq, respectively, for smooth potentials and he found by using a limit in which 
the smooth potential reduces to a hard core potential that for hard spheres 
«4 = 0(|i| 5 ) and kq = 0(\t\ 7 ), while the odd cumulants were still zero. An 
alternative approach was followed by De Schepper et al.{3| consisting of 
directly evaluating the Van Hove self-correlation function for short times for 
hard spheres based on pseudo-Liouville operators (which replace the usual 
ones for smooth potentials). De Schepper et al.|25j obtained n n = 0(|i| n+1 ) 
for all even n > 2, with corrections of 0(|t| n+2 ). 

The 0(t 2n ) result for smooth potentials is the more remarkable in that 
a naive estimate of the short time behavior of K n based on its connection 
with the moments would predict a behavior as the n-th order in t. Hence all 
terms from 0(t n ) up to 0{t 2n ~~ 1 ), 1 should vanish. The question we address 
here is whether this is indeed general for smooth potentials. 

Everywhere in this paper "up to 0(t a y means "up to and including 0(t a )'" . 
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In fact in this paper we shall show that under quite general conditions, 
the main one being that in the initial ensemble, the velocities have a multi- 
variate Gaussian distribution |2J] and are statistically independent of the 
positions, one can prove a very general theorem to be presented in the 
next part. This theorem implies, among other things, that for a many 
particle system in which the (interparticle as well as external) forces are 
independent of the velocities, the nth order cumulants of the displacement 
are indeed 0(t 2n ) for n > 2 and 0(t n ) for n < 2. This confirms the above 
stated expectation, but extends it to a restricted class of non-equilibrium 
situations as well, as we will show. 



1.2 Class of physical systems 

In this part we will restrict ourselves to the following class of systems (the 
Theorem in Part |2] of the paper concerns a more general class). Consider 
N point particles in three dimensions 2 whose positions and velocities are 
denoted by the three-dimensional vectors rj and Vj, respectively, and whose 
masses are mj, where i = 1...N. The time evolution of the particles is 
governed by 

i-i = Vj (1.3) 
v< = FiCr* t)/mi, (1.4) 

where Fj is the force acting on particle i, which is supposed to be a smooth 
function of the time t and of r , which is the collection of all positions rj. 
Likewise, will denote the collection of all velocities Vj. 

We consider an ensemble of such systems. In the ensemble, the initial 
probability distribution V(r N , v N ) of the positions and the velocities of the 
particles is such that the velocities each have a Gaussian distribution and 
are statistically independent of the positions, i.e., 3 



N 



r(r N ,v N )=f(r N )U 



A™A 3/2 ( i fl i 

exp - -Pimi\Vi - Uj 



2ir 1 V 2 



(1.5) 



where Uj is the average of the initial velocity Vj and Pi are positive "inverse 
temperature" -like variables. In Eq. H1.5JI . the f(v N ) denotes a general prob- 
ability distribution function of the initial positions of the particles. While 

2 Other dimensionalities work just as well. 

3 The expression in Eq. 11.51 is not the most general Gaussian distribution for v . 
While the Theorem in Part [5] of the paper allows any multivariate Gaussian distribution 
of the velocities, we have in fact not found any physical applications for that yet. 
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the distribution function in Eqs. (|1.5j) can describe both equilibrium and 
non-equilibrium situations, it always shares with the equilibrium distribu- 
tion the Gaussian dependence on the velocities, which is in fact crucial for 
the Theorem below to hold. 

We stress that the forces Fj in Eq. (|1.4jl are independent of the veloci- 
ties, but may depend on the positions and on the time t in any way as 
long as they are smooth. An example of smooth forces would be infinitely 
differentiable forces Fj (r N , t) , but also Lennard- Jones- like forces are allowed 
provided the distribution of the positions f(r N ) assigns a vanishing prob- 
ability to the particles to be at zero distance of one another, which is the 
singular point of the Lennard-Jones-like potential at which it is not smooth. 

Canonical equilibrium for a single or multi-component fluid is just a 
special case of these systems. In that case, one has Uj = 0, 0i = (3, Fj = 
—dU/dvi and f(r N ) oc exp[—f3U(r N )], where U{r^) is the potential energy 
function of the system. 

The probability distribution functions of the form Eq. (|1.5(l in which 
each particle has its own mean velocity Uj and "inverse temperature" Pi, may 
seem at first of a mathematical generality which has little physical relevance. 
Note, however, that this is a convenient way to describe mixtures of any 
arbitrary number of components. In such a mixture, the mean velocities 
and/or temperatures of the different components could be selected physically 
e.g. by having two vessels with different substances at different temperatures 
with a divider which is opened at t = 0, or by means of a laser (or perhaps 
even a neutron beam) applied at t = 0, tuned to a resonance of one of 
the components only, which would give the particles of that component a 
nonzero average momentum as well as a different initial "temperature" due 
to the recoil energy [l{J. To what extent such techniques could ensure the 
Gaussianity of the velocity distribution is a technical matter that we will 
not go into here. 

1.3 Definition of moments and cumulants 

In this paper, when we speak of a random variable, we mean a function of 
the positions and the velocities and possibly the time t. All physical 
quantities are variables of this kind. For such variables A(r , v N , t), the 
average with respect to the distribution function V(r N , v N ) will be denoted 
by 

(A) = J dr N dv N P(r N ,v N )A(r N ,v N ,t). (1.6) 
Throughout the paper, A's will denote general variables. 
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We define Arj(i) = (Axj(i), Ayj(i), Azj(i)) as the displacement of parti- 
cle i in time t, where i = 1 . . . N and the time dependence is suppressed for 
brevity. The nth moment 4 of Axj(i) is the average of its nth power, and is 
denoted as 

Mn = (A<(i)), (1.7) 

where the dependence on i and t on the left-hand side (lhs) has been sup- 
pressed. In the following we will also suppress the t dependence in Axj. The 
cumulants of Axj, 4 denoted by equal to the moments \x n with certain 

factorizations of them subtracted and thus sensitive to the correlations of 
Axj. Their precise definition is via the cumulant generating function |2^. 

oo 

n! 

n=l 



$(fc) = log^exp[ifeAxi]) =X)^( ifc ) n - ( U 
from which the K n follow as 



Hn = 



d(ik) n 



(1.9) 



fc=0 



An alternative notation for the cumulants, which is more analogous to 
E q. (II. 7|) and which is especially convenient in the case of several variables, 
is 24] 

«A*J n] ))^ n . (1.10) 

We stress here that the superscript [n] is not a power. Furthermore, we will 
follow the convention that the superscript will not be denoted if n = 1, i.e., 
{{Axi)) = ((Ax\ )) (which is also equal to k\ = (Axi)). 

The cumulant generating function in Eq. (|1.8|) can be expressed in terms 
of the moments fx n since $(fc) = log(exp[i/cAa;j]) = log[l + ^^i(i^) n A i n/^!]- 
The relations between the cumulants and moments can then be found by 
Taylor expanding the logarithm around 1 and using a multinomial expansion 
for the powers of the second term inside the logarithm. This gives, somewhat 
formally, 



ill 



OO OO / / fl]\Pl 

{pe>0} 1=1 e=i ll ' 



4 Although the moments fi n and cumulants k„ are really the "moments and cumulants 
of the probability distribution function of Axi(t)" we will refer to them simply as the 
"moments and cumulants of Axi(t)." 
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For instance, for the first few k„., Eg. becomes «i = /Ui, «2 = A*2 — Mi 

and k 3 = fj, 3 - + 2^0,13. 

In general, K n = /i n =t factored terms, where the factored terms contain 
all ways of partitioning the moment fj, n into a product of lower moments \ig 
such that all I values (taking into account their "frequencies of occurrence" 
Pi) add up to n. 

Apart from the cumulants n n of the single variable Axj , we will also need 
the general definition of cumulants which applies to the displacement of the 
same particle in different spatial directions as well as the displacements of 
different particles. Like for the cumulants of a single particle's displacement, 
the expressions for the general cumulants are moments with factored forms 
subtracted, e.g. 

((Asi;As 2 » = (AxiAx 2 ) - (Axi)(Ax 2 ) (1.12) 
((AsijAtftjAzi)) = (AxiAyiAzi) - (Axi)(AyiAzi) - (AxiA yi )(Azi) 

-(AxjAsiXAjft) + 2(Ax 1 )(Ay 1 )(Az 1 ) (1.13) 
<(a4 2 W 2] >> = (Ax 2 Ay 2 ) - (Ax 2 ) (Ay 2 ) - 2(Ax!A yi ) 2 
-2(Ax!)(Ax!Ay 2 ) - 2(Ax 2 A yi )(A yi ) 
+8(Ax 1 )(Ax 1 Ay 1 )(Ay 1 ) - 6(Ax!) 2 (A yi ) 2 . (1.14) 

Here, the semi-colons are inserted to avoid ambiguity and to make it explicit 
that the expressions inside the double angular brackets are not to be mul- 
tiplied. This is different from the notation for the multivariate cumulants 



used by Van Kampen |2J] who denotes the above cumulants as ((AxiAx 2 )), 
((Axi AyiA,zi)) and ((Axf Ay 2 )), respectively, which would suggest that the 
expressions inside the double brackets are products of powers, which they 
are not. 

The general definition of the multivariate cumulants will be given in 
Eq. (|2.9|) in Sec. 12.11 (where we will also show how the various prefactors are 



determined) and can also be found in refs. [23| and |2 
1.4 The Theorem in physical terms 

For the class of physical systems described above, the general and somewhat 
formal Theorem in Part [21 of the paper states that, in a less abstract physical 
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notation 



(( Ax[ ni * ] ; Ay [ ^ ] ; Az [ ^ ] Ax l ? 2x] ; Ay [ ^ ] ; Az [ ^ ] ;•••)) 

= (c {niv} t" + 0(t^) ifn<2 
\c {n „ )} t 2 " + 0(t 2 "+ 1 ) if n>2, 

where are non-negative numbers (some may be zero), % = 1 . . . N, rj = 
x, y or z and 

N 

n = Y^ n "? ( L16 ) 

i=l rj=x,y,z 

is the order of the cumulant on the lhs of Eq. (|1.15|) . 

Note that the cumulants n n that occur in the expansion of the Van Hove 
self-correlation function are special cases of these nth order cumulants, i.e., 
we can write 



c n t zn + C(r n+1 ) for n > 2. 



Section [2.41 contains the expressions for the c n and the C{„.}. 

We stress once more that on the basis of the connection between the 



cumulants and the moments (jl.llj) . only an 0{t n ) scaling of the nth order 
cumulants is to be expected, so that this is a nontrivial theorem. The result 
in Eq. (|1.17f) generalizes the results of Schofieldj^] and Sears 0] who found 
that for an equilibrium liquid, K2 = 0{t 2 ) while K4 = 0(t 8 ) and kq = 0(t 12 ). 
The results in Eqs. Q1.15JI and (jl.l7j) hold, however, also out of equilibrium 
as long as the initial distribution of the particle velocities is Gaussian and 
statistically independent of the particles' positions. 

1.5 Applications 

We will now discuss a number of physical applications of the Theorem in 
Eqs. (|1.15|l - (|1.17jl just presented. 

1.5.1 The equilibrium Van Hove self-correlation function 

a) In incoherent neutron scattering on an equilibrium fluid, one essentially 
measures the equilibrium Van Hove self-correlation function G s (r,t) (with 
r = |r— r'|), which is the Fourier inverse of the incoherent scattering function 
K,(M)IJ B 0, ED, EH, El- If the length of the wave vector k in F s (k,t) 
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is small, then according to Eq. Q1.2J1 F s (k,t) ~ exp[— K 2 k 2 /2], i.e., nearly 
Gaussian, and so its inverse Fourier transform G s (r,t) is also approximately 
Gaussian. Corrections to this Gaussian behavior can be found by resumming 
F s (k,t) = exp^2 n K n (ik) n /n\ in Eq. JOJ) to the form 



F s (k,t) 



e -\n 2 k 2 



1+ ]T ^ 



n=A 
n even 



(1.18) 



using that in equilibrium odd cumulants are zero. For n > 4, the coefficients 
b n are given bv[17|| 



e n \^m" 

{Pi>0;i even} £=4 



;i.i9) 



£ even 



For example, 64 = K4/4!, &6 = kq/61, bs = (ks + 35k|)/8!. Fourier inverting 
the resummed form of F s (k,t) leads to 



G s (r,t) 



1 + 



00 



11-20) 



Here -ff n is the rath Hermite polynomial and the dimensionless w = r /^J2k 2 . 
We note that the series in Eq. (|1.2U|) appears to have a fairly rapid 
convergencejy]. Taking just the first few terms would give 



G s (r,t) 



\J2lXK 2 



1 



6!8k: 



;i.2i) 



The Theorem in Eq. 1)1.17)1 provides a justification for the expansion in 
Eq. (|1.20() of the self part of the Van Hove function G s for short times t in 
the following way. As Eqs. 1)1. 19|) and 1)1.20)1 show, the cumulants /-c n >4 give, 
via the b n , corrections to a Gaussian behavior of G s (r,t). The Gaussian 
factor e~ w suggests that typical values of w are 0(1) in Eq. 1)1.20)) . so also 
H n (w) =C(1). Its prefactor in Eq. ()1.20)) is, however, t-dependent through 
b n /(2K 2 ) n l 2 . Given the relation between b n and n n in Eq. ()1.19)) it is easy 
to see that they scale similarly, i.e., if the conditions of the Theorem are 
satisfied so that K n>2 = 0(t 2n ) then also b n>2 = 0{t 2n ). Since k 2 = 0(t 2 ), 
we obtain 

, K /0 = 0(t n ). (1.22) 
(2At 2 ) n/2 
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This means that the series in Eq. (jl.2Uf) is well-behaved for small times 
t, in that each next term is smaller than the previous one, and that by 
truncating the series 5 one obtains for small t approximations which can be 
systematically improved by taking more terms into account. Note that in 
contrast if n n had been 0(t n ), each term in the series in Eq. ()1.2U)) would 
have been of the same order. 

b) An expansion of a similar form as Eq. 1)1.20)1 was found by Rahman p| 
for G s (r,t), and by Nijboer and Rahman[5( for F s (k,t). Their expressions 
are in terms of the so-called non- Gaussian parameters a n . These non- 
Gaussian parameters have recently also found applications in the context 
of supercooled liquids and glasses, where the y h ave been proposed as a kind 



of order parameter for the glass transition |l ll . 112 . Ilq and as measures of "dy- 
namical heterogeneities" in supercooled liquids and glasses (Isl. Efll. 21. E3|. 
We note that while supercooled liquids and glasses are not in true equilib- 
rium (for that would be the solid phase), they do have a Gaussian velocity 
distribution 'inherited' from the fluid phase. 

Given the present interest in these non-Gaussian parameters a n , we will 
now compare them with the cumulants n n . The non-Gaussian parameters 
a n are defined in terms of the distance AR = \J Ax 2 (t) + Ay 2 (t) + Az 2 (t) 
traveled by any particle in time t in a three dimensional fluid, as [J] 

(AR 2n ) 

Un = (Ai?2)«(2n + l)!!/3« (L23) 

For isotropic fluids, (AR 2n ) = (2n + l){Ax 2n (t)) = (2ra + V)H2n, so that 
Eq. ()1.23j) can then be written as 

_ ^ 2w -(2n-l)!!^ 
n ~ (2n-l)!!/iS ■ 1 ' 

We now see that even though both the a n and the cumulants n n >2 are, by 
construction, zero for Gaussian distributed variables, in Eq. (|1.24j) the a n 
are 2n-th moments \i2n with only the most factored term, fi^, subtracted, 
while the cumulants K2 n in Eq. (|1.11|) have all possible factored terms sub- 
tracted. Using Eq. 1)1.24(1 and the inverse of the relation between K n and fj, n 
in Eq. 1)1.11)1 (which can be found by using the generating functions), it is 

5 When truncating, terms that are kept have to include not just the leading order in t, 
but at least all terms up to 0(t k ) if terms up to Hk are kept. In numerical approaches it 
is possible to retain the full terms. 
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possible to express the a n in terms of 



a,. 



e n 

0<p e <n 1=1 



pt 



11-25) 



According to this formal relation, the first few a n are given by 



«2 



«3 



1 K4 
3 K 2 



K i 



+ 



15 «| 



(1.26a) 
(1.26b) 



The Theorem in Eq. (|1.17l) says that for small times t, k 2 = 0(t 2 ) and 
Kn>2 = 0(t 2n ), so that k 2 „/k2 = 0(t 2n ). According to Eq. (fOKJl . all a„ 
have a contribution from k^/k 2 , so that all a n are of 0(i 4 ), in contrast to 
the cumulants K 2n , which are of increasing order in t with increasing n. In 
fact, using Eq. (|1.25|) and the Theorem, one can derive straightforwardly 
that the dominant term for small t in Eq. H1.25j) is the one with p\ = n — 2, 
P2 = 1 and pi>2 = 0, which leads to 



ft, 



n(n — 1) 



«2 



1.27) 



plus a correction of C(t 6 ). Thus, for small t, 03 is approximately three times 
ot2, «4 six times a 2 etc. Such approximate relations are indeed borne out 
by Rahman's data on a 2 , 03 and 0C4 for small t0, 0|. 

In terms of the a n , the expansion of G s in Eq. (|1.21|) becomes, with the 
help of Eq. (fPfil . 



G,(r,t) 



1 i 2,012 u ( \ _i 15 (" 3 ~ 
1+ 4!4^ H + 6!8 



3«2 



-i?e(^) + 



1.28) 



Although formally equivalent to Eq. (|1,21|) . for small t, Eq. Q1.28JI is some- 
what less convenient, because one cannot see right away that the last term 
(involving Hq) is 0(t e ) rather than 0(i 4 ), as one may naively suspect from 
a 2 = 0(t A ) and 03 = 0(i 4 ). Thus, a cancellation between the leading orders 
of 03 and 3a 2 has to take place, which, especially using a numerical evalu- 
ation of the a n , may be hard to obtain. On the other hand, in Eq. i|1.21|) . 
this cancellation is automatic. 

We note that, somewhat surprisingly, the data of Rahman also show that 
the relation between a n and a 2 in Eq. (|1.27|) is still approximately satisfied 
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2) hard sphere-like 




-0.1 ps ~1 ps -10 ps 

t 



Figure 1: Sketch of the behavior of the non-Gaussian parameter ai = 
K4/(3k|) f° r regular liquids in equilibrium in different time regimes (based 
on fig. 7 in ref. |4| and fig. 4.11 in ref. [8|), with our physical interpretation 
for each regime. 

for larger times 0]. This seems even true for hard spheres [2(|. Thus the 
higher order non-Gaussian parameters a n >2 are apparently dominated by 
«2 for larger times just as they are for smaller times t. This dominance of a.2 
makes it hard to extract from these higher order non-Gaussian parameters 
any information that was not already contained in 02- The cumulants K, n , 
or perhaps the b n , may contain additional information about correlations in 
(supercooled) fluids in a more accessible form (compared to the a n ), and 
may therefore be a more suitable choice to investigate such correlations for 
all times t. 

c) Returning to the series in Eq. (|1.20() . although it may be well-behaved 
for small enough i, it is not known up to what time this remains so. Since 
b n = 0(t 2n ) only to leading order in t, there is a point in time after which 
the 0(t 2n ) term is no longer a good approximation, and this may be related 
to the point at which the series in Eq. (|1.20j) is no longer guaranteed to be 
useful. Rahman investigated 02 (as well as 03 and 04) numerically for a 
model of liquid argon (temperature 94 K, density 1.4 • 10 3 kg/m 3 )|4|. (See 
also ref. 0] for a broader overview.) Figure 1 shows a sketch of the non- 
Gaussian parameter ct2 as a function of i, based on fig. 7 in ref. |J{ and 
fig. 4.11 in ref. 0. One sees that 02 is a very flat function near t = 0, which 
persists only up to roughly t ~ 0.1-0.2 ps. At that point the curve shoots 
up rapidly, leading to a large "hump", which lasts up to about 10 ps (or 
perhaps somewhat below that), after which it starts to decrease to zero. 
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Although somewhat outside the scope of this paper, we would like to 
give a possible interpretation of the numerical results sketched in Figure 1 
for moderate and high densities. 1 ) The flat behavior of near t = corre- 
sponds to the 0(t 4 ) behavior as given by the Theorem in Eq. Q1.17[) . 2) Be- 
cause hard spheres can be seen as a limit of a smooth interparticle potential 
in which the steepness goes to infinity 0], and the potential used by Rahman 
is rather steep, the shoot-up phenomenon at ~ .1 p s may well be related 
to the hard-spheres result of De Schepper et al. 25] that K n = 0(\t\ n+1 ), 
or a2 = 0(\t\ 5 /t 4 ) = 0(\t\), as follows. The steep but smooth potential of 
Rahman will resemble a hard sphere fluid on time scales t s on which a col- 
lision has been completed [23]. Thus at t = t s the scaling 0(t A ) for 02 ought 
to go over to 0(|i|), or in fact, because of the duration t s of the 'collision', 
to 0(\t — t s \), which would require the kind of sharp increase observed by 
Rahman at ~ 0.1 ps. 3) While a persistence of non-Gaussianity occurs 
already in dilute systems [2^] . this is due to a different mechanism than in 
denser systems, where it comes about because the particle is trapped in a 
"cage" formed by its neighboring particles, with which it has repeated and 
correlated collisions. This effect is dominant at high densities 4) The 
decay of 02 to zero indicates that the motion becomes Gaussian and pre- 
sumably sets in (for dense systems) when the particle manages to escape 
its cage. After escaping, it finds itself in a new cage environment consisting 
largely of particles with which it has not interacted before. This motion from 
cage to cage is called cage diffusion (2^. 5) From a central limit theorem 
argument using that successively visited cages after many cage escapes have 
little correlation with each other, one would then expect Gaussian (and pre- 
sumably but not necessarily diffusive) behavior. The precise mechanism of 
the long-time behavior at low densities falls outside the scope of this paper. 

Also in simulations of a supercooled argon- like mixture[l9|. «2 plotted 
as a function of time shows a flat curve for short times and a sharp increase 
around 0.1 ps, while «3(i) shows similar behavior. The interesting part 
from the perspective of supercooled liquids and glasses, however, is in how 
far that increase continues and on what time scale and how 02 decays back 
to zero, which takes a very long time for supercooled liquids and is related 
to the time scale at which particles escape their cages. But the Theorem in 
Eq. (|1.15j) -Eq. (|1,17|) has nothing to say about ct2 on that time scale. 

1.5.2 Local equilibrium systems 

As a first example of a non-equilibrium system to which the theorem may 
apply, we consider a fluid not too far from equilibrium, such that it has ini- 
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tially roughly an equilibrium distribution except that the temperature, fluid 
velocity and density are spatially dependent, i.e., it is in local equilibrium. 
The class of initial distributions in Eq. (jl.5|) does not seem to be of that 
form, and indeed, if Pi and Uj are allowed to vary with rj, then the proof 
as given in the next part of this paper, runs into difficulties. Nonetheless, 
one can construct distributions of the form ()1.5j) which physically describe 
precisely the local equilibrium situation. Imagine dividing the physical vol- 
ume V into M cells of equal size and assigning to each cell a a temperature 
Pa, a fluid velocity u a and a density n a . The particles of the system are 
divided up as well, putting N a = n a V/M particles in each cell, such that 
particles 1 through JVi are in cell 1, N± + 1 through N\ + N2 are in cell 2 
etc. This can be accomplished by choosing /(r ) in Eq. (|1.5|) such that the 
chance for these particles to be outside their cell is zero. Next, we set all 
the Pi and Uj of the particles in cell a equal to p a and u a . If the cells are 
big enough so that fluctuations in the number of particles can be neglected, 
this situation describes local equilibrium just as well as spatially dependent 
P(ri) and u(rj) can, and for this constructed local equilibrium, the Theorem 
in Eq. (jTTTHj) holds. 

1.5.3 Out-of-equilibrium phenomena on very short time scales 

In this section, we consider somewhat more general out-of-equilibrium sys- 
tems, namely those that start at t = with a Gaussian velocity distribution. 
The distribution of the positions of the particles is allowed to be very dif- 
ferent from that in equilibrium, however. 

The cumulants whose short time scaling was obtained here also occur 
naturally in the Green's function theory, which was developed for far-from- 
equilibrium phenomena on the picosecond and nanometer scales [lil. E3. Iljj. 



Ilfitll7l|. Considering for example a mixture of two components, one can write 
the density of component A (A = 1 or 2) at position x and for simplicity 
here in one dimension as: ^3] 

n x (x,t)= dx'Gx(x,x',t)n x (x',0), (1.29) 



where the non-equilibrium Green's function G\(x,x' ,t) is the probability 
that a particle of component A is at position x at time t given that it was at 
position x' at time zero. By an expansion detailed in a future publication [l^| 
G\ can be written similarly as G s in Eq. (|1.2()j) . as 



G x (x,x',t) 



1 + V- ^—HAw) 

£- 3 (2 K2 )«/2 > 



(1.30) 
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(1)1 



Here the dimensionless w = (x — x' — k{)/^2k2 and 

(L;!I) 

{p e >o} e=3 

Some examples are: 63 = K3/3!, 64 = K4/4!, 65 = and &6 = 

(k6 + 10k|)/6!. These equations are non-equilibrium generalizations of the 
equations for the equilibrium Van Hove self-correlation function in (|1,2U|) 
and Eqs. Q1.19JI . respectively. However, different from the case of the Van 
Hove self- function, in the Green's function theory, the K n , and thus the b n 
through Eq. (|1.31|) . depend on x' because in that theory the single particle 
i of component A is required to have been at the position x' at time zero. 
This requirement can be imposed by multiplying the probability distribu- 
tion function in Eq. (jl.5ft by 5{xi — x') (times a proper normalization). The 
resulting distribution describes the subensemble of the original ensemble for 
which particle i is at x' at time zero. Note also that it is still of the same 
form as Eq. Q1.5JI . with /(r ) — > f'(r N ) = 5{xi — x')f(r N ), where x' is just 
a parameter. Thus, as long as we restrict ourselves to systems with Gaus- 
sian initial velocity distributions, the Theorem in Eq. H1.17j) still applies. 
Note that the fix ) may describe any non-equilibrium distribution of the 
positions of the particles. 

The Gaussian factor in Eq. (|1.3U|) suggests as before that typical values 
of w are 0(1), so H n (w) = 0(1). Its prefactor in Eq. (fOHl) is b n /(2K 2 ) n/2 . 
Using the Theorem in Eq. (|1.17|) that n n = 0(t 2n ) and the relation between 
b n and n n in Eq. (|1,31|) it is easy to see that also b n = 0(t 2n ). Since 
K2 = 0(t 2 ), we obtain 

, 6 " l0 = 0{t n ). (1.32) 

This means that, for these systems, the series for the non- equilibrium Green's 
functions G\ in Eq. (jl.30f) are well-behaved for small times t (just as the 
equilibrium Van Hove self- correlation function was) and that by truncating 
the series one obtains for small t approximations which can be systematically 
improved by taking more terms into account. This statement is in fact not re- 
stricted to the one-dimensional case of the Green's functions discussed here; 
the three-dimensional version, for which the b n become tensors and which 
will be presented elsewhere |l7||. exhibits — when the general formulation in 
Eq. (|1.15|) is used — the same scaling with t of the subsequent terms in the 
series for the Green's functions. 
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We note that in non-equilibrium situations, «2 = ^4/(3^2) may have 
a similar behavior as sketched in Figure 1 for the equilibrium 012. Al- 
though a proper numerical test is yet to be performed, this expectation 
is roughly consistent with numerical results of the Green's function for heat 
transport [3, [H, OH. In that case the contribution of the non-Gaussian 
corrections in Eq. 1)1.30(1 (involving the Hermite polynomials) were most sig- 
nificant on the sub-picosecond time scale, whereas extrapolation indicated 
that hydro dynamic-like results may occur for times possibly as short as 
2 PS0. 

In view of the possible application to nano-technology, it is important 
to understand the behavior on very short time scales and the related small 
length scales over which a particle typically moves at such time scales. The 
Green's function theory can potentially describe a system on all time scales. 
The current Theorem assures that this theory can at least consistently de- 
scribe the very short time scales, by showing that the expansion of the 
Green's function is well-behaved. For practical applications, and to know 
how short the time scales must be for the Theorem to apply, it is still neces- 
sary to determine the coefficient of the (D(t^ n ) of K n: i.e. the c n in Eq. 1)1.17(1 
(cf. Sec. 12.4(1 . This will require a numerical evaluation of the moments of 
derivatives of the forces, which we plan to do in the future. The behavior of 
K n at longer time scales, will also be investigated in the future, with special 
attention to the question whether the Green's functions give hydrodynamic 
behavior for long times 17 ] . 



2 Mathematical Part 

In this part, we present (Sec. 12.2(1 and prove (Sec. 12.3(1 the Theorem in its 
most general, mathematical form. 

In the mathematical formulation of the Theorem, as well as in its proof, 
it is convenient to adopt a different notation than the one introduced in 
Part Since from a mathematical point of view it does not matter whether 
different degrees of freedom are associated with different spatial directions 
or with different particles, the mathematical notation will treat all these 
degrees of freedom on the same footing. This can be achieved by calling the 
x component of particle one the first degree of freedom, its y component the 
second, its z component the third, and then the x component of particle two 
the fourth degree of freedom, etc. The total number of degrees of freedom 
is then M = 3N. 

One can achieve the goal of describing all degrees of freedom in the same 
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way by associating with each degree of freedom a (generalized) position n 
and a (generalized) velocity Vi, and rewriting the three dimensional real 
positions and velocities Vj by using the mapping 

























w 



(2.1) 

and similar for the velocities and the average velocities, i.e., 

vi \v 2 | , v 2 U 5 | , ••• (2-2) 




vi -> \v 2 , v 2 -»■ U5 , ... (2.3) 



respectively. 

Applying this mapping on the equations of motion Eqs. 1)1.3(1 and 1)1.4(1 . 
one finds that the generalized positions ri and velocities Vi (collectively de- 
noted by and v^ respectively) satisfy 

h = Vi (2.4) 
Vi = ai{r N ,t) (2.5) 

where the a,i(r~^,t) are the accelerations which follow from the mapping 

U(r",t) , U^.t) , ...(2.6) 

Wr«,t)J m2 W^)7 

We note that a possible interpretation of this formulation is that of J\f 
"particles" moving in just one spatial dimension. This convenient picture, 
in which each and V{ may be seen as the one-dimensional position and 
velocity of a "particle" i respectively, will be used below, even though the 
system is not really one-dimensional. Cumulants that in the physical for- 
mulation of the previous part involved both different spatial directions as 
well as different particles (as e.g. on the lhs of Eq. 1)1.15(1 ). become in this 
picture cumulants involving just different "particles" (degrees of freedom), 
and these kinds of cumulants are therefore the quantities of interest here. 

We remark that the proof of the general Theorem that we will give below 
is a "physicists' proof", meaning that the proof does not claim to have com- 
plete mathematical rigor but has every appearance of being correct, perhaps 
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under mild and reasonable additional conditions such as the existence of the 
moments and cumulants. 

First, however, we need to introduce the general definition of multivariate 
cumulants [23|, [2J], needed in the Theorem and its proof, and some of their 
properties. 



2.1 Preliminaries: the general definition and properties of 
cumulants 

Here we will give the general definition of cumulants which applies to any 
number of variables, i.e., of the multivariate cumulants ijj 2]. Thereto we 



introduce a set of (general) random variables A q with q = 1 . . . Q, for which 
we assume that a probability distribution function exists. The cumulant 
generating function of the A q is defined as 



. . . kg) =log(^expJ2 ik <i A i)- ( 2 - 7 ) 

g=l 

Here ( ) denotes an average with the probability distribution function of the 
A q . Note that if the A q are expressed in terms of yet other random variables 
(later, in the Theorem, the and v^) the average may equivalently be 
taken with their probability distribution function. With the help of this 
generating function, the cumulants can be expressed similarly as in Eq. (|1.9|) 
by 




«A^ ; ...;4jd) 

{fej=0 



(2. 



Note that {k q } = is a short-hand notation for k\ = 0, hi = 0, . . . , fcg = 0. 
Furthermore, when one of the n q is equal to 1, we will omit the corresponding 
superscript [1], e.g., ((^i;^)) = ((A^;A^]). We stress once more that our 
notation is different from Van Kampen's|2J], who writes ((A™ 1 A 1 ?, 2 . . . Aq)). 
In particular, the square brackets in the superscripts are intended to show 
that they are not powers of the A q , while the semi-colons in the expressions 
inside the double angular brackets indicate that one should not interpret 
them as products, but as defined by Eq. ()2.8|) . Note also that the expression 

fn 1 

A q q only has a meaning inside a cumulant. 

The cumulants can be expressed in terms of moments analogously to 
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-ni\---n Q \ Y, Em" 1 ! II 



1 / (A?...A$) \pw 

^ y^-r, ) 11 | I 4!-..^ ! J ' 

T,{i} i qP{t}=n q 

(2.9) 

where {£} = {£%, . . . ,£q} denotes a set of Q nonnegative l q values, q = 
1 . . . Q and p^ gives the frequency of occurrence of that set. In Eq. (|2.9|) 
the sum is over the frequencies p^y of all possible sets {£} (and thus all 

possible moments (A^ 1 . . . Aq)) for which Y^=i ^qP{e} = n q- This is nothing 
else than factorizing the expression A™ 1 ...A n Q Q in all possible ways. Hence, 

the cumulants ((A^ l] ;. . . ; A [ q q] )) = (A™ 1 . . . A n Q Q )± factored terms that are 

products of the moments (A^ 1 . . . Aq ) such that the l q values for fixed q add 
up (taking into account their frequencies Pyey) to n q . For some examples of 
multivariate cumulants, see Eqs. (|1.12l) - (|1.14(l . 

We furthermore define the order n of a cumulant as the sum n = 
Y^q=i n q- Note that the maximum number of moments that are multiplied 
in any term in Eq. (|2.9j) is (taking into account the frequencies) equal to the 
order of the cumulant. 

We will now give four properties of the cumulants that follow from its 
definition in Eqs. (|2.7|) and (|2.8|) and that can be used to manipulate ex- 
pressions inside the double angular brackets of the cumulants. 

(a) The first property is that, as with averages, constants may be taken in 
front of cumulants: 

(((CA 1 )^;...))=C^((A [ ^ ] ;...)) (2.10) 

(b) The second property concerns cumulants where the same quantity occurs 
more than once: 

((4" l] ;4 n2] ;--.» = «4 ni+n2] ;•••)) (2-11) 

(c) The third is a multinomial expansion: 

(((A 1 + ...+^)W ; ...» = £ -^-((4^;...;^;...)). 

{k q }>0 llg=l fe 9- 
Y.%k q =n 

(2.12) 
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(d) Finally, cumulants have theproperty that they are shift invariant if their 
orders are larger than onej^, Hjj- E.g., if C\ is a constant, 

«(A 1 + C7 1 )W ; ...» = ((AW;...». (2.13) 
Only if the order of the cumulant is one, does a shift have any effect: 

((A 1 + C 1 )) = ((A 1 )) + C 1 . (2.14) 

2.2 The Theorem 

Theorem. For a classical dynamical system of N degrees of freedom de- 
scribed by (generalized) coordinates r, and velocities V{ (i = 1...M), col- 
lectively denoted by and respectively, for which 

a) the time evolution is given by the equations of motion 

h = Vi (2.15) 
Vl = ai{r M ,t) (2.16) 

with velocity-independent, smooth accelerations aj, and 

b) the initial ensemble is described by a probability distribution in which 
the velocities are multivariate Gaussian 24] and independent of the 
coordinates r%, i.e., with a distribution function of the form 

M 



p(r^,iT) = /(r Ar ) x /det [E/(2vr)] exp 



- ^ijiVi-Ui^Vj-Uj) 



2 



(2-17) 

where — is a positive, symmetric J\f x J\f matrix with constant ele- 
ments Hjj, itj is the average velocity corresponding to the ith degree of 
freedom, and f(r^) is the probability distribution of the coordinates. 

When these conditions are satisfied, the cumulants 6 

« tlH} = ((Arl Bl ';Ar^;...;Ar^» (2.18) 

of the displacements 

Ar, = n{t) - n(0) (2.19) 
6 These cumulants are denned in Sec. 12. II with Q — > TV, q — + i and A q — > Ari 
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(whose t dependence has been suppressed) satisfy for sufficiently short initial 
times t 

= {c {ni} t n + 0(t n+1 ) ifn<2 
Hni] \c {ni} ^ + 0(^ +1 ) ifn>2, 
where n is the order of the cumulant given by 



(2.20) 



n 



n\ + n 2 + . . . + nj^f. 



(2.21) 



The coefficients C{ n .j will be given later in Eq. (|2.44|) in Sec. 12.41 



2.3 Proof of the Theorem 
Strategy based on Gaussian velocities 

The Theorem formulated in Eqs. H2.15tJ2~2*Tj) in Sec. 12.21 will be proved in 
this section, although we will defer the details of the proof of a required 
auxiliary theorem to the Appendix for greater clarity of the procedure. 

To obtain the initial, short time behavior of the moments and cumulants 
of the displacement, Ar^ may be Taylor-MacLaurin expanded around t = 
as 



Arc 



P 



dt"f t=o 7! ~~ 4rf 



7=1 



dPi- 



t=o 7! 



(2.22) 



where we used (|2.15j) . Because of the equations of motion (|2.15|) and (|2.16|) , 

the d^Vi/df 1 , viewed as functions ol 
itly of t, are recursively related by 



the d J Vi/dt~i ', viewed as functions of and (at time t) as well as explic 



d^ 



dt-r+ x 



E 



d /d?Vj 



Uj dr- 



dt"f 



d fd^Vi 



dv, V dV 



+ 



d /d^Vi 



dt V df 



(2.23) 



To show that the cumulant K{ ni } 



0{t ) for n > 2, it is of course 
possible to work out this cumulant straightforwardly using Eqs. (|2.9j) . (|2.22|) . 
(|2.23f) and ()2.17|) . in that order. Such a procedure was essentially followed by 

SchofieldH for K n = K{ n fi,o,...} = ((^ r i"')) for n < 6 and SearsQ for n < 8 for 
equilibrium fluids. In their expressions many cancellations occurred before 
K n could be seen to be, for these cases, of the 2n-th order in t instead of the 
n-th order in t. These cancellations seemed to happen as a consequence of 
equilibrium properties. However, by carrying out the same straightforward 
procedure for the more general class of non-equilibrium initial conditions 
in Eq. ()2.17|) . we have found that, while odd moments are no longer zero, 
still K n = 0(t 2n ) for n = 3, 4, 5 and 6. These results for K3, K4, and 
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k,q naturally led us to propose the Theorem. Because the straightforward 
calculations for this non-equilibrium case are very lengthy, they will not be 
presented here. In any case this procedure is not very suited to determine the 
order in t of ft{ ni } for general rii, because with increasing order n = Yli=i n i 
an increasing number of terms have to be combined (taking together equal 
powers of t from the various products of moments) before they can be shown 
to be zero. 

Our strategy for proving that nth order cumulants K{ ni \ = 0(t 2n ) for all 
n > 2, will be to exploit the Gaussian distribution of the velocities as much as 
possible. But we can only hope to use the Gaussian nature of the velocities 
if we succeed in bringing out explicitly the dependence of the coefficients of 
the power series in t of the cumulants on the velocities. This dependence 
has so far only been given implicitly — the cumulants are related to the 
moments by Eq. Q2.9|) . the moments contain [Arj] , Arj is expanded in the 
time t in Eq. (|2.22jl , and the coefficients in that expansion are the derivatives 
of the velocity i>j, which can be found by using Eq. H2.23j) recursively. To 
make this more explicit, the first step of the proof will be to expand the 
cumulants as a powers series in the time t and the second step will be to 
express this series more explicitly in the velocities. In the third and last 
step we will then use the properties of Gaussian distributed variables, i.e., 
the velocities, to complete the proof of the Theorem. 

It turns out that the properties of Gaussian distributed variables that 
we will require in the third step of the proof are formulated for independent 
Gaussian variables whose mean is zero, while in Eq. (|2.17|) the velocities are 
Gaussian but do not have zero mean, nor are they independent (because 
3jj ^ for i ^ j in general). For this reason, it is convenient to intro- 
duce already at this point new velocity variables whose mean is zero (cf. 
Eq. tfFm ): 

Vi = Vi- (vi) =Vi- m. (2.24) 

The Vi are generalizations of the peculiar velocities used in kinetic theory 
and will be referred to in the general case treated in this paper as peculiar 
velocities as well. For the same reason, it is convenient to get rid of the 
statistical dependence of the initial velocities in Eq. (|2.17|) by bringing the 
(positive) matrix E to its diagonal form £i5jj. This can be accomplished 
by an orthogonal transformation. We will assume here that this orthogonal 
transformation has been performed so that E is diagonal, and will show at 
the end of the proof that the form Eq. (|2.20j) of the Theorem is invariant 
under such a transformation (while the coefficients C{ n .} do change). 

Substituting Eq. (|2.24|) into the time expansion in ()2.22j) and into the 
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probability distribution Eq. (|2.17|) (using that E is now diagonal), gives 
the expansion of Ar^ and the probability distribution V{r^ ,V^) in their 
peculiar velocity form: 



A? 



u; 



' + ^ 7! dt-y- 1 
7=1 



t=o 



(2.25) 



V(r^,V M ) = /(r^nv^expl-i&lf] (2.26) 

i=l 

respectively. In order to make future expressions less complicated we intro- 
duce for the coefficients in Eq. H2.25j) the notation (7 = 1,2,...) 

X i Jr M ,V M ) = ^^^ . (2.27) 

7V 1 7! din- 1 t=o y ' 

Thus Eq. ((2~23|) becomes 

00 

Ar i = u l t + Y,X ll {r M ,V M )t\ (2.28) 
7=1 

We remark here that in the double indices of Xi~(r", V^), the first index 
always pertains to a particle number and will be denoted by the roman 
letter i, while the second pertains to an order in t and will be denoted by 
the Greek letter 7. Below, we will drop the explicit dependence of the Xj 7 
on r M and V N . 

We will now start the actual proof of the Theorem for general n. 



First step: Expanding the cumulants in the time t 

The infinite number of terms in the time expansion Ar^ in Eq. (|2.28j) 
means that to get the power expansion in t of the cumulants H{ ni } = 

((Arj™ 1 ^ ; . . . ; Arjy- )) occurring on the lhs of Eq. (|2.2U|) . we would use 
Eq. lESJ) with 

Q -»■ M, q^i, A q ^ An. (2.29) 

and combine the infinite number of terms coming from Eq. (|2.28l) . Obviously 
if we are interested in the cumulants up to 0(t 2n ~ 1 ), where n = Y^iLi n i 
(cf. Eq. ()2.21[0 . we should not have to retain all these terms in Eq. ()2.28() . 
but only those up to 0(i 2n_1 ). As a matter of fact, we need even less terms, 
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namely only terms up to 0{t n ) in Eq. (|2.28|) . as the following reasoning 
shows. The cumulants occurring on the lhs of Eq. (|2.20f) are given in terms 
of the moments (Arf 1 Ar e 2 2 . . . Arty) = (WiLi Ar f> b y E Q- UM- Taking 
the terms up to the n-th order in t in Eq. (|2.28|) . i.e., writing Arj = Uit + 
X)"=i X^t 1 +(D(t n+1 ), it is straightforward to show that the moments satisfy 

i=l i=l 7=1 

so that, for given pm, 

N M n 

n ( n = n ( n m+e v)*) + ° (tn+E m Pw ] - 

{1} i=l {1} i=l 7=1 

(2.31) 

Products of this kind occur in the definition of the cumulants on the rhs of 
Eq. (|2.9[) . and are summed over with the restriction that (with q = i 

here) ^{f} £iP{e} = ^i- Since also Yli=x n i = n > Eq. <|2.31j) becomes 

N AT n 

n(n A ^r } = u(u^ t +E x ^y i ) PW+o ^- (2 - 32) 

{Z} i=l {Z} i=l 7=1 

As this holds for each such expression on the rhs of Eq. Q2.9JI (with Eq. I|2.29|) ). 
we have for its lhs 

K{ni} = ((Ar[-l;...;Ar^)) 



( Ul t + J2 M ; • • • ; (tw* + E x -^ 7 ) [rwl )) 



7=1 7=1 

+0(t 2n ). (2.33) 

We remark that the first term on the rhs of this equation gives all powers 
t n up to t 2n ~ 1 , which we are interested in, as well as some of the powers 
of t higher than 2n, which we are not interested in. The second term, i.e., 
0(t 2n ) only contains powers of t of 2n and higher. So with Eq. (|2.33f) we 
have established that for the lower powers of t only n terms in the time 
expansion of Ar^ are needed, but we have not separated the powers of t 
lower and higher than 2n completely yet. 

For that purpose we continue from Eq. (|2.33jl and use first the shift 
invariance property of cumulants explained at the end of Sec. 12.11 to obtain 
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from Eq. (l2~33j) 



(X;jri7* 7 ) NI ;-;(EW) 

7=1 7=1 
AT 

j=i 



(2.34) 



Furthermore, using the multinomial expansion in Eq. (|2.12j) . we can write 



XV) N ;...;(EW) 



[njv] V 



7=1 



E 

{n^ 7 >0} 



A' 



n 



7=1 



i n™=i 



<(xft ll] ; Xg 12] ; . . . ; £"=i «*r7 



(2.35) 



Note that the summation indices n- n arise from the multinomial expan- 
sion, where i denotes a particle index and 7 runs from 1 to n. For ft{ rai } in 
Eq. (|2.34|) . we need this quantity only explicitly up to 0(i 2n_1 ), so powers of 
t higher than 2n — 1 may be discarded (i.e., combined with the 0(t 2n ) term) 
and only powers lower than 2n need to be kept. Hence, in the exponent on 
the rhs of Eq. (|2.35|) . we only need terms with Yli=i ^7=1 n «77 < 2n. Since 



also n = YA=iT^=i n i-n we obtain TZ=i n i-yl < 2 T, J i=iT^=i 



A' 



n 



or Yli=i n ii- > Ya=i S"=2 n «7(7 ~~ 2). Combining this condition with 
Eqs. (1231) and (|23ffl) . and using X a = V t (cf. Eq. STTTt ). we find that 
the expansion of an nth order cumulant in time t up to 0{t 2n ~ l ) is given by 



-A' 



A' 



n 



1 n"=i 



n 



17- 



E 

{ni 7 >0} 

E^inii>E^iE^ = 2"i7(7-2) 



A' 



.2n\ 



(2.36) 



n<5 nj0 +o(t 2 

i=l j^i 

We note that in Eq. (j2.36f) only cumulants appear, instead of moments, 
and, more importantly, that powers of t lower and higher than 2n are easily 
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identified, something that in the straightforward moment-based approach 
mentioned above only happens after a lengthy calculation. 



Second step: Writing cumulants in terms of peculiar velocities 

The dependence of Kt ni \ in Eq. (|2.36|) on the peculiar velocities V™ follows 
from the dependence of the X,^ on the V . 

According to Eq. (|2.27|) . the coefficients JQ 7 can be determined using 

X^,V^) = X iy (r^,V^,t)\ t=Q , (2.37) 

where we have defined 



**r(r^>*) = ^^ (2-38) 



taken at time t. Below, we will suppress the explicit dependence of X^ on 
r~", and t. The X^ satisfy the recursion relation, from Eqs. (|2.23() : 



X, 



17+1 



7 + 1 



N 

E 

3=1 



( i \r\®X%-y , j\f ,dXi-y 

( u j + v j)-^r + a 3( rJ ^~dV~ 



(2.39) 



For 7 = 2, Eqs. (|2~3%j) . (l2~21) and (|2~TH|) show that X i2 = (dVi/dt)/2 = 
t), which is independent of and is thus a polynomial in the pecu- 



liar velocities of total degree zero. Using the recursion relation Eq. (|2,39|) 
one sees that if X^ is a polynomial in of total degree 7 — 2 then on the 
rhs of Eq. ()2.39|) . the term (uj + V^dX^/drj is a polynomial of total degree 
7 — 1, while ajdXiy/dVj has a total degree of 7 — 3, and dX^/dt has a total 
degree 7 — 2. The highest total power of in Xj 7 +i is thus 7 — 1. 

In other words, for 7 > 2, is a polynomial in the peculiar velocities 
V N of total de gree 7 — 2, with coefficients that can depend on the positions 
of the particles. From Eq. ()2.37|) we then see that also X^ = X^t = 0) is 
a polynomial in the peculiar velocities of total degree 7 — 2. Therefore 
we can write for the dependence of X^ on V 

Xi, = E E h {Pj} (r^)Vr---V^ (2.40) 
p=°T,jPj=p 

It turns out that in the third step of the proof we will not need the precise 
and explicit forms of the b^ p .y(r^), but only that the total degree of the 
polynomial X^ is 7 — 2. 
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Third step: Using the Gaussian nature of velocities 

Given the polynomial nature of the function of V N in Eq. (|2~4TH) . wc 

can give the following interpretation to the second condition under the sum- 
mation sign in Eq. (|2.36|) . i.e., Y^Li n n > Si=i Y^=2 n h^ ~ 2 )- Eacn ex ~ 
pression of the form x\^ ] in the cumulant ((v} nn] ; . . . ; V$ Ml] ; X l ™ 12] ; . . . ; X [ £\ 
in Eq. (|2.36|) (although devoid of meaning outside of cumulant brackets) 
signifies that the highest power of X^ occuring inside the aver- 

ages in the expression for the cumulant in terms of moments on the right 
hand side of Eq. (|2.9j) . Because Xj 7 is a polynomial in F-™ of total de- 
gree (7 — 2), this highest power of X^ is a polynomial in of total 

In- 1 

degree 71^(7 — 2). Then, all such expressions -X" i7 ' 7 together in the cu- 
mulant ((V} nil] ; . . . ; Vp m] ;X [ £ 12] ; . . . ; X [ ^ n] )) in Eq. generate poly- 
nomials of at most degree X^7=2 n n(7 ~~ 2 ) inside the averages in the ex- 
pression for the cumulant in terms of moments on the right hand side 
of Eq. Q2.9JI . On the other hand the expressions in the cumulant 
{(V} nil] ;. . . ; Vp m] ;X [ £ 12] ; . . . ; X [ ^ n] )) in Eq. together generate poly- 
nomials of at most degree YliLi n ii- Thus, the condition YhLi n n > 
Yli=i Ylj=2 n n (7 — 2) in Eq. (|2.36f) indicates that the total degree in 

generated by the V- in the cumulant {{v} nil] ; . . . ; Vp m] ; X [ ™ 12] ; . . . ; X [ J% n] )) 
is larger than the total degree generated by the X^. 

The crucial point is now that the auxiliary Theorem lAl in the Appendix 
can be applied to cumulants of this form. This theorem is most conve- 
niently expressed in terms of general random variables P q , which were used 
there as well to denote general random variables, but which are now poly- 
nomial functions of the peculiar velocities V . Theorem ^states that if P q 
(q = 1 . . . Q) are polynomials of total degree d q in the independent, zero- 
mean Gaussian distributed V^, and n q and pi are nonnegative integers, 

then ((V^ ; . . . ; V$" ] ; Pj ni] ; . . . ; P [ q q] » = if Pi > £« =1 n q d q , except 

when all n q >i = and one pi = 2, in which case it becomes just ((V^ }}. 

In order to apply Theorem ^ to each term in the summation on the 
rhs of Eq. tCT . we need to rewrite ((vf* 1 ,. . . ; V$ m] ; X l £ 12] ; . . . ; X [ ^ n] )) 

as ((V} pi] ;. . . ; V^ 1 ; p{ m] ; . . . ; Pq q] )). This can be achieved by a mapping 
of single to double indices, i.e., by setting 

Q->M(n-l) (2.41a) 
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(Pl,ni,di)->-(Xi2,ni2,0) (P n ,n n ,d n ) pf 22 ,n 22 ,0) 

(-P 2 ,ra 2 ,d 2 ) -> (Xi3,ni3, 1) (P n+ i,n n+ i,d n+ i) -> (A 23 ,n 23 , 1) 

(P 3 ,n3,d 3 )^(Xi4,ni4,2) : 



(P n _i,n n _i,d n _i) -> (Xi n ,ni n ,n - 2) (Pq,uq, <Iq) -> (Xj^f n ,nj^ n ,n - 2) 

(2.41b) 

and 

Pi -> "a (2.41c) 

Then Y^=i n q^q ls seen to be equal to S™ =1 ^7(7 — 2) = d. As 

the restriction on the summation in Eq. (|2.36f) shows, for all terms on the 
rhs of Eq. (|2.36|) (except the 0(t 2n ) of course), Y^a=i n n ^ s larger than 
this d. Theorem tells us that the cumulant occurring in each term of 
these terms is then zero except when nj 2 = = ■ ■ ■ m n = for all i 
and only one nn = 2 and rij^n = 0. This exception means, since also 
2_4=i X^7=i n n = n -> that n = rin = 2. So the only possible nonzero term 
in the sum in Eq. (|2.36|) occurs for n = 2. Furthermore, for n = 1, the last 
term Y^i=i u i*<W Yij^i ^0 in Eq. (|2.36|) is also left. 

Thus, using Theorem lAl we have shown that each term in Eq. (|2.36|) is 
zero separately except for n = 1 and n = 2, so that 

(0(t n ) if n < 2 , 
|ni} \o(i 2 ™) ifn>2 ^ ' 

remains on the rhs of Eq. (|2.36|) . Given that Kr ni \ can be expanded as a 
power series in t, Eq. 1)2.42)1 coincides with the formulation Eq. 1)2. 20j) of the 
Theorem, if we define C{ n .} as the coefficients of the t n and t 2n , respectively. 
This is therefore now proved (with the proviso that Theorem A is proved 
in the Appendix), but up to this point only for diagonal matrices E in the 



distribution 1)2.17)1 . 

To obtain the same result for non-diagonal matrices E in Eq. ()2.17)) . we 
apply a transformation S (an orthogonal J\f x J\f matrix with elements Sij) 
such that E' = S • E • S T is diagonal. In this transformation, r[ = Sj=i <% r j> 
v[ = Ylj=iSij v j and also Ar. = Y^=i^ij^ r j- Since =! is diagonal, the 
cumulants of the primed displacements Ar^ satisfy the Theorem. The cu- 
mulants of the original displacements Ar^ = Ylj=i Sji^ r j can be expressed 
in terms of the primed ones using the multinomial expansion 1)2.12)1 . with 
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the result 



E m) _ _ nj^l 
tj.i 1 I • • • n-i »rl ^— ' 



nn! • • • nij\fi ^ nj\f\\ ■ ■ ■ nj\fj\f\ 

{nij} {nUji 

x(((S 11 Ar[)^;...;(S m Ar' x )^; 

(5 12 Ari) [n2l] ; • • • ; (5^ 2 ArV) [ "^ ] ; . . . ; (S^ArV)^ 1 )) 

(2.43) 

Since in each term on the right hand side the n *i = Si=i n « = n 

(cf. Eq. (|2.21j) ). and each term contains the cumulants of primed displace- 
ments, which were already shown to satisfy the Theorem, each term scales 
as t 2n if n > 2 and as t n if n < 2, and therefore so does the sum. This proves 
that Eq. Q2.2UJI of the Theorem holds for arbitrary (positive symmetric) ma- 
trices E. □ 

2.4 Expression for the coefficients in the Theorem 

We will now show how one can determine the coefficients of the t 2n term of 
K {ni}i i-e. c {ni} in Eq. (|2.20[). 

For n = 1 and n = 2, it is straightforward to show that 24] 

C{i,o,o,o,...} = = «i (2.44a) 

c {2 ,o,o,o,..} = (^i 2 ) = E^ln ( 2 - 44 b) 
c { i,i,o,o,..} = = [- _1 ]i2- (2.44c) 

(and similarly for c {0j i A o,o,-} and c {o,2,o,o,o,-}' c {i,o,i,o,o,-}' c {o,i,l,o,o,-}' 
etc.). 

To find C{ n .} for n > 2 one can use a similar calculation of K( ni \ as 
used above, but taking one additional term Xi n +it n+1 in the time expan- 
sion of Arj in Eq. ()2.33[) into account. Performing then the same kind of 
manipulations as in the proof above, one arrives for n > 2 at 



{n i7 >0} 

^ -y=l n i-y =n i 



ri7=i "'- 



77 f 



X ((A n , ...,A 1 +1 ,A 21 , ...,A N +1 )). (2.44dj 
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where n was defined in Eq. I|2.21jl . 

Note that for the case of a cumulant of one displacement, e.g., ((Ar™)) 
in Eq. (|1.17|) . one has 

Cn = C{... i0 ,n,0,...} ( 2 - 4 5) 

where the n on the rhs is at the ith position. 

Using Eq. (I2.44dj) . we can give some examples of c n for n = 3 and n = 4: 



c 3 = 3((V^;X t4 ))+6((V t ;X i2 ;X i3 )) + ((X^}} (2.46) 
c 4 = 4((V^;X l5 }} + 6((v} 2] ;X\l ] }} + 12((v} 2] ;X t2 ;X t4 )} 

+12((V i ;X\ 2 };X i3 )) + ((X§)). (2.47) 

while an example of the coefficients C{ n .} for cumulants involving different 
degrees of freedom, e.g. ((An; Ar 2 ; Ar 3 )) is, from Eq. (|2.44d|) . 

c { i,i,1A0...} = ((Vv, V 2 ;X U )) + ((Vr,X 22 ;X 33 )) + ((F i; X 23 ; X 32 )) 

+ ((Vf, X 24 ; X 31 }) + ((X 12 ; V 2 ; X 33 )) + ((X 12 ;X 22] X 32 }) 
+ ((X 12 ; X 23 ;V 3 )) + ((X 13 ; F 2 ; X 32 )} 

+ ((X 13 ; X 22 ; V 3 )) + ((X u ; V 2 ; V 3 )). (2.48) 

We note that although the X^ are useful to derive these expressions for 
C| n .}, to evaluate them in practice requires additional work. One would first 
need to write the JQ 7 as X^ir^ ', , t = 0) using Eq. (|2.37|) . and work out 
the recursion relation ()2.39j) . Furthermore, the cumulants would have to be 
written in terms of averages using Eq. (|2.9|) . The values of these averages 
will depend on the system, i.e., the accelerations Oj, and their evaluation 
will in general require a numerical approach, which we will explore in future 
work. 



Conclusions 

In this paper, we have proved a mathematical theorem on the correlations 
of the initial time displacements of particles (or in general of coordinates 
associated with degrees of freedom) for a class of dynamical systems whose 
main restriction is that the initial distribution of the velocity-variables is a 
(multivariate) Gaussian, independent of the position- variables. 

Among the physical applications of this Theorem is the result that a well- 
known short-time expansion of the Van Hove self-correlation is well-behaved: 
each subsequent term in the expansion is smaller than the previous one for 
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small enough t, something which had been suspected but not established 
before. It has also been shown on the basis of the Theorem that in the studies 
of undercooled liquids and glasses, using cumulants instead of the usual non- 
Gaussian parameters may give more physical information. Furthermore, 
it was shown that the expansion used in the Green's functions theory is 
also well-behaved if the velocity distributions are initially Gaussian, so that 
this theory, which can describe non-equilibrium mass transport processes on 
short time and length scales, has now been given a firmer basis. 

Acknowledgments 

The authors would like to thank Prof. H. van Beijeren and Prof. F. Bonetto 
for carefully reading the manuscript and for their helpful suggestions. This 
work was supported by the Office of Basic Energy Sciences of the US De- 
partment of Energy under grant number DE-FG-02-88-ER13847. 

Appendix 

In the third step of the proof of the main Theorem, we needed the following 
auxiliary theorem to prove the main Theorem. 

Theorem A. Let V* be a set of N statistically independent, zero-mean 
Gaussian variables, collectively denoted by V^, and let P q (q = 1 . . . Q) be 
a set of Q polynomials in of total degree d q . Let pi (i = 1 . . . M) and n q 
(q = 1 . . . Q) be nonnegative integer numbers. Then if 




(A.l) 



i=i 



q=l 



and at least one n q / 0, the following cumulant vanishes: 



//T/[pi]. . T/bw]. pb 

\\ v l > • ■ ■ > V M > "\ 



;P^ ] » = 0. 



(A.2) 



while if all n q are zero, one has 




-;^»=x;«^»^ii^o. 



(A.3) 



i=l j^i 
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Proof of Theorem{^ Since we have not found a proof of this theorem in the 
literature, we will give it here, but before we can prove Theorem ^ we need 
four lemmas and the definition of a (^-modified average. This definition 
will serve, in conjunction with Lemma ^ below, to construct a convenient 
generating function (which takes the form of a ^-modified cumulant) for 

the quantities ((v} Pl] ;. . . ; ; P[ ni] ; . . . ; P [ q q] }) that occur in Theorem \B 
Working with this generating function will be more convenient than try- 
ing to calculate each ((V^ l] ;. . . ; V$ x \ Pj" l] ; . . . ; P [ q q] )) individually. After 
this generating function has been introduced, Lemma [21 and Lemma |3] are 
presented and proved, which use the Gaussian nature of the V{ to give prop- 
erties of averages and ^-modified averages of powers of the V% , which allow 
one to prove a polynomial property (Lemma|3J) of the generating function of 
the quantities ((v} Pl] ;. . . ; V$" ] ; p\ ni] ; . . . ; pjg q] }}. This polynomial property 
will be the central ingredient to complete the proof. 

First we remark that below, will always denote the same set of 
AT statistically independent, zero-mean Gaussian distributed variables Vi of 
Theorem ^1 

The definition of ^-modified averages and cumulants is: 

Definition. The ^-modified average of a general variable P (which below 
will always be a polynomial in V^) is defined as 



where the 9i (i = 1 . . . jV) are real numbers. 

Similarly, ^-modified moments are generally defined as the ^-modified 
averages of powers of general variables (i.e., functions of ', and in the 
main text also of r^), and ^-modified cumulants are defined as having 
the same relation to ^-modified moments as normal cumulants have to 
normal moments, ^-modified cumulants are therefore also given through 
the ^-modified cumulant generating function as: 

Q 

The generating function nature of the ^-modified cumulants follows from: 

Lemma 1. (relation between cumulants and 9^ -modified cumulants) 7 . The 
#-^-modified cumulants of a set of variables Pi , ... Pq are related to the 

7 For J\f = Q = m = 1 this lemma coincides with the last exercise of section XVI. 3 in 
Van Kampen's book|24l| 



«P^ ; . . . ; P^ )) e , = J] 7^ ^ ( ex P E ^ 
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normal cumulants by 



((P|" lJ ;...;P^V=E- E 

Pl=0 PM=0 



n 

.i=l 



x((V^...;V^Pt l] -,...;P$ Q] )) (A.6) 
if at least one n q 7^ 0, while it is zero otherwise. 

Proof. For case in which at least one n q is nonzero, we start with the rhs of 
Eq. HA,6|) and use the expression for the cumulants in Eq. H2.8|) : 



Pi=0 PjV=0 



6>? ! 



n „,i 



i=l 



//T/[Pll. . T/[PAr]. pM. . pK?l\ 

\\ "l i • • • j v Af > M ' " • • ' Q / 



00 co M aPt pa>i Q on, 

p 1= PA r=0i=l ^ v * ; q=l y q/ 



M Q 

lOg ( [ ^ \k[Vi + ^ ikgPg 



i=l 



9=1 



(A.7) 



{M={*0=° 



We recognize the Taylor series, i.e., that Y^=o([~^i\ Pi / P^-) d Pi f(k[)/dk 
/(— iOi) to write this as 



k'=0 



~M aPi 

TT -i- 

11 Pi! 

pi=0 PaT=0 Li=l 



E-E 



//T/[Pil. . T/[PAr] . p[«l] . . p[ n Ql\ 
\\ "l i • • • J *V > M ' ' ' ' ' Q / 



Q f)n q N Q 

II 5^ io«(«p[E^+E^ 



1=1 



9=1 



(A.8) 



{fc,}=0 



Using definition (|A.4|) . we obtain 
•JV 



E-E 

Pl=0 PjV=0 



n 

i=l Pl \ 



//T/[pi]. .i/biv].pM. . pI n Q]\ 

\\ "l ' * * * ' > -M ' * * " ' Q / 



n 



d n 1 



l a(ifcq 



<3 -, N 

g=l 1=1 



(A.9) 



{fc 9 }=0 
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Since we are considering the case that at least one n q is nonzero, the contri- 
bution from the /cq-independent \ X)t=i ^l^Yi) vanishes when d/d(ik q ) acts 
on it. The remainder on the rhs of Eq. (|A.9|) is by definition ()A.5|) equal to 
the lhs of Eq. (lA~6|) . 



The case n% = ■ ■ ■ = uq = is trivial, for then ((-P™ 1 . . . Pq Q ))qm = since 
a zeroth cumulant is always zero. 

□ 



A consequence of this lemma, i.e., of Eq. (|A.6|) . is that repeated deriva- 
tives with respect to 8{ of ((P-j ni '; . . . ; Pq Q ~^))qM taken at 6>j = generate the 

((^ 1 ];...;y^];P^;...;it Q] ))- 



The second lemma concerns the average of a product of powers of the Vt 
and a single polynomial in . 

Lemma 2. (Averages of powers ofVi times a polynomial in V ) Let P be a 
polynomial function of total degree d in '. Given a set of M nonnegative 
integers Sj, collectively denoted by s-^ , and a set of N "parities" <5j with 
each 5i zero or one, one can write 

( II V^p) = P s * ({ Sj }) H(2 Si + 2* - W(V?y* (A.10) 

i=l i=l 

where P s *({sj}) is a polynomial in the Sj of total degree s*, which satisfies 

d-^r N s- 

s* < 2 ■ (A ' U) 
Proof. We start by writing out the polynomial P as 

P= E h P >}fl V ? 1 - (A.12) 

{**} «=i 

Z? =1 Pi<d 

Of these terms, only those with the same "parity" as the 5i (i.e., if <5j = 0, pi 
is even and if 5i = 1, Pi is odd) contribute to the average (n^=i V^ St+&1 P) , 
due to the even nature of the distribution of the V{. Thus, one can write 
Pi = 2s\ + 5i and obtain 

fiv^p)= £ 6 {2s:+M n^ 2(si+s ^ i) )- ( a - 13 ) 

i=l ' {s'J i=l 
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Using (V? s ) = (2s - ) s |23j, we obtain 

i=l ' {s'J i=l 

= E Wall {2Si + 26i _ m 

K> »=1 V 1 ' 

Af 

x Y[(2 Si + 25 t -l)\\(V?) s * (A.14) 
i=i 

This is of the form stated in Eq. (jA.lOj) . where 

pam) = e w*> n n ( 2s ^ + + ^ - i)(^ 2 ) <+di 

E^i(2<+5 t )<rf 

(A.15) 

This shows that P s * is a polynomial in , since it depends on the Sj only 
through finitely many factors (2sj + 2s'- + 25j — 1). In fact, for each term in 
Eq. (|A.15|) . the number of such factors is Y^iLi s 'i- Because of the restriction 
on the sum over {s'j} in Eq. (|A.15|) . this number YliLi s i 1S ^ ess than or 
equal to (d— Yli=x $i)/2 for each term. The total degree s* of the polynomial 
P s * ({sj}), which is the maximum of this number over all the terms, therefore 
also satisfies s* < {d - Y^L 1 Si)/2. □ 

The third lemma concerns a polynomial property of the 0^-modified 
averages of polynomials P in V . 

Lemma 3. (polynomial property of the 6^ -modified average of a polynomial 
in ') Let P be a polynomial of total degree dmV M . Then the 6 M - 
modified average of P is a polynomial in 9^ of at most the same total 
degree d. 

Proof, a) Consider first the case that P is a polynomial in the Vi of total 
degree d and of definite "parity" for each i, i.e., that it either changes sign 
or remains unchanged when Vi is replaced by — Vi. Define the numbers Si 
such that 5i = if P is even in Vi and 5i = 1 if it is odd. By the definition 
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of the ^-modified average in Eq. I|A.4|1 . expanding the exponent, and using 
that only even powers of Vi have a non-zero average, one gets 



e 2 

oo oo 

xE-E 

si=0 sj\f=0 



(28 i + 5 i )\ 



Y[V? Si+Si P) (A. 16) 



Using Eq. (|A~Tol) of Lemma 01 and that (2s;-l + 25M\/(2s,+5M = l/(2 s * Si l), 
we obtain 



{P) 6 M = e 2 



' Af 



i=l 



{ Si >0} i=l 

ls^Af /)2/,/2\ 



2 s ' s A 



{ Sl >0} 



Note that the polynomial P s * in Sj inside the {si} summation can be gen- 
erated from an expression that does not have this polynomial by applying 
operators 9 2 A8/89 2 j ) } e.g., 



E.n« 

{ Sl >0} i=i % - 



P>Q2 2^ 11 



9 2 

J { Sl >0}i=l 



9 1 v-AT 



3 oel 



£2 



Ef =1 ef<^> 



(A.18) 



and in general 

e ^({,}>n [ ^ 2r 



Ps*({0]^})e^^). (A.19) 



{s;>0} i=l 

Combining this with Eq. (|A. 17|) gives 



(A.20) 

In this expression, the differential operators (8/89?) "bring down" factors of 
#| from the exponent to its right. That exponent itself is then canceled by 
the exponent to the left of the P s * operator, so that only a polynomial in 9^ 
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is left. The total degree of this polynomial is twice the maximum number 
of factors brought down by the operators (twice because the squares of the 
6j are brought down). This maximum number is simply the total degree 
of P s *, i.e., s*. Counting finally also the powers of 9i of the first product 
in Eq. (|A.2fl|) . we see that {P)$m is a polynomial in 9^ of total degree 
d! = 2s* + YaLi^- But Lemma H in particular Eq. (IA~TT1) . says that 
s* < (d — Y^Li £j)/2, so that for a polynomial P of definite parity, its total 
degree d' satisfies d 1 < d. 

b) To show that the same is true for a general polynomial, note that 
any polynomial P can always be written as a sum of polynomials of definite 
parity. Since the ^-modified average of that sum is the sum of the 9 M - 
modified average of each term, and each term is a polynomial in 9^ of total 
degree d! < d, (P)gN is also a polynomial in 9^ of total degree d! < d. □ 

The next and final lemma we need before we can prove Theorem 
concerns ^-modified cumulants of several polynomials in ' . 

Lemma 4. (9^ -modified cumulants of polynomials in ) Let P q (q = 
1 . . . Q) be a set of polynomials of total degree d q in the V . Then the 

^-modified cumulant ((Pj ni] ; . . . ; Pq q] )) 0jv is a polynomial in 0™ of a total 

degree of at most d = Y2q=i n gdq- 

Proof. Because the ^-modified cumulants are defined formally in precisely 
the same way as normal cumulants [i.e., Eq. ()A.5|) vs. Eq. (|2.8|) ]. the relation 
Eq. I|2.9|) between cumulants and moments applies to the 0^-modified cu- 
mulants and moments as well. According to that relation, the #-^-modified 
cumulant ((pj m] ; . . . ; P^ q] ))qM can be expressed as a sum of terms each of 

which contains a product of the moments (P^ 1 ■ ■ ■ P^)gM raised to the power 
P{i} (where the product is over all possible sets {£} = {li, . . . , Iq} with the 
restrictions stated in Eq. (|2.9|0 . Because P q is required by the conditions of 
the lemma to be a polynomial in V of total degree d q , P q q is a polynomial 
m V N of total de gree £ q d q , and the product P^ 1 ■ ■ ■ Pq is a polynomial in 
V N of total de gree Y^=i^qdq- According to Lemma |3 each ^-modified 
moment (P^ 1 ■ ■ ■ Pq Q )qm is then a polynomial in 9^ of a total degree of at 
most ^2q = i^qd q . Its p^yth power in Eq. ()2.9() is then a polynomial in 9^ 
of total degree P{i}Y^=\^qdq ( a t most), and the product (over {£}) that 
occurs on the rhs of Eq. (|2,9|) is of total degree d = Yl<{e.}P{t}Y^=i^qdq 
(at most). Since p^ is summed over in Eq. (|2.9|) with the restriction that 
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Yl{e}P{e}^q = n <?> the total degree d of the expressions n^j^i 1 ' ' ' ^Q 3 )^^ 
can be rewritten as d = Ylq=i n qdq (at most). Each term in the sum over 
on the rhs of Eq. (|2.9j) is therefore (at most) of this total degree d, so 

that the full expression, i.e. ((p]™ 1 '; . . . ; Pq Q ^))qm, is also a polynomial in 
9^ of a total degree of at most d. □ 

Conclusion of the proof of Theorem XM 

We can now finish the proof of Theorem To consider ((vj Pl ^ ; . . . ; vj^^ ; 
p] ni ^ ; . . . ; Pq Q ^)) as on the lhs of Eq. (|A.2|) in the case that at least one n q ^ 
0, we look at its generating function, i.e., ((p[ ni ^ ; . . . ; Pq Q ^))oM. According to 
Eq. (|A.6|) of Lemma^ this generating function admits a power series in 
At the same time, according to Lemma 0] (whose proof required Lemmas |2] 

and|2J), ((-P]" 1 '; . . . ; PQ Q ^))e^ ls given by a polynomial in 9^ of a total degree 
of at most d = Ylq=i n qd q , so its power series in 9^ in Eq. (|A.6|) terminates 
(at the latest) after YliLiPi = Since these lemmas should be valid for 
all 9^, each coefficient of the terms with Y2i=iPi > ^ must be zero, i.e., 
using Eq. JOJ) of Lemma □ ((v} pi] ;. . . ; ; p[ ni] ; . . . ; P [ q q1 )) = 0, which 
is Eq. HA.2|) of Theorem lAl Finally, for the case n q = for all q = 1 . . . Q, 
Eq. HA.3|) of Theorem lAl follows simply from the zero-mean Gaussian nature 
of the Vi and their statistical independence. □ 
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